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Abstract 

We derive a general crack propagation law for slow brittle cracking, in two 
and three dimensions, using symmetry, gauge invariance, and gradient expan- 
sions. Our derivation provides explicit justification for the "principle of local 
symmetry," which has been used extensively to describe two dimensional crack 
growth, but goes beyond that principle to describe three dimensional crack phe- 
nomena as well. We also find that there are new materials properties needed to 
describe the growth of general cracks in three dimensions, besides the fracture 
toughness and elastic constants previously used to describe cracking. 



I. Introduction 

There are many aspects of the problem of crack growth that have received a 
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lot of attention recently. For instance, there has been much interest in dynamic 
fracture-*^ and the accompanying crack bifurcation^'^ and other instabilities.^ 
The transition between failure due to percolation of a network of many small 
cracks and failure due to a single dominating crack has also been explored,^ as 
well as the transition between brittle and ductile cracking.^ Pattern formation 
in multiple cracking^ has also been of interest. In light of all the interest in 
these rather complex phenomena of fracture, it is somewhat surprising to find 
that little is known about the growth laws for even slow-growing, simple (not 
multiple) three dimensional cracks, though there has been some work done on 
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calculating the paths of cracks in two" and three dimensions, and many 
measurements and calculations of the crack velocity for simple two-dimensional 
geometries. "'^'^^"'^^ 

The problem of finding a growth law for cracks would seem to be of fun- 
damental interest; so, in this paper, we apply the standard tools of theoretical 
physics — gradient expansions, symmetry, and gauge symmetry — to find the most 
general possible growth law for a three dimensional crack growing slowly in a 
homogeneous, isotropic medium. Since it is possible to make precise numeri- 
cal computations of the elastic fields for arbitrary three-dimensional geometries, 
with today's computers, in a matter of hours, we do not consider the related 
problem of finding the stress state of the material containing the crack, but 
consider it to be completely known. We also compare the crack growth law we 
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derive here to previously derived and measured properties of cracks in two and 
three dimensions. In a second paper, we will discuss the detailed behavior of 
cracks growing under this law, using linear stability analysis as well as numerical 
simulations; in a third paper we will discuss experiments designed to measure 
material-dependent parameters appearing in the crack growth law. 

II. Simplifications 

We begin by simplifying the problem of crack propagation using length and 
time scale considerations. First, we smooth our crack problem over the length 
scale Is which characterizes the size of inhomogeneities and anisotropies in the 
material containing the crack. (For example, in a glass, is is a few atomic 
sizes; in a polycrystal, it is the grain size; in concrete, it is the size or distance 
between the pebbles it contains.) Although for a single crystal, is is as large as 
the body containing the crack, for many situations of practical relevance, ig is 
much smaller than the size of the body. In those cases, we can smooth the crack 
problem over is without losing much information, making the crack a smooth 
surface, and the material containing the crack continuous, homogeneous, and 
isotropic. 

A second length scale in crack propagation problems arises because every 
material has some stress above which it fails to have linear elastic properties (e. g. 
begin to flow, with plastic or viscous behavior; emit dislocations; break bonds; or 
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have a martensitic transformation). For some materials, this stress is very low, 
and there is no linear elastic regime at all. For others, linear elasticity is valid 
except very near the crack tip, where the stress is much higher than in the bulk 
of the body. For these materials, there is a length scale £ni which characterizes 
the size of the non-linear process zone around the crack tip; £ni can range from 
a few angstroms in glass to tens of centimeters in concrete. In this analysis, 
we consider only materials for which £ni is small compared to the length of the 
crack and the size of the body, so that the bulk of the material can be considered 
linear elastic. This work, then, describes materials usually considered linear and 
brittle, as well as materials exhibiting viscoelasticity, plasticity, and martensitic 
transformation toughening, as long as the length scale for these behaviors is 
sufficiently small. In principle, the non-linear properties of such materials could 
be included in a later version of this work to extend its applicability to smaller 
length scales. 

A third length scale relevant to crack propagation is associated with the 
degree of translational invariance along the crack front. For many crack systems 
studied in the past, every plane perpendicular to the crack front is equivalent, 
which means that the problems can be considered two dimensional. On the 
other hand, most practical crack problems are not two dimensional, but instead 
have crack front curvature or stresses which vary along the crack front. If this 
is the case, then there is a length scale, which we call the dimensional crossover 
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length Icici above which the problem is three dimensional; f^c is either one of 
the geometric lengths associated with the crack geometry (such as the radius of 
curvature of the crack front), or is associated with the stress gradient: ~ 
For this work, we assume that f^c is large, though in general not as large as the 
size of the body containing the crack, and we expand in powers of quantities 
which are inversely proportional to (,cLc (i- e. gradients). 

Finally, we also simplify the crack propagation problem by considering 
cracks which are growing slowly enough that inertial and relativistic (relative 
to the sound, not light, velocity) effects are unimportant. Some of these effects 
could be included in future work, but the present analysis suffices for cracks 
which arrest after growing a certain distance, such as when a wedge is driven 
into a crack; cracks which grow at a constant speed, such as under constant 
displacement loading; and cracks which may eventually speed up, but which are 
currently growing slowly, as in the cases of fatigue cracks, sub-critical cracking, 
and the first stages of growth under constant force loading. 

III. Relevant variables 

Although the knowledge of length scales from the previous section has sim- 
plified our problem to a nearly two dimensional, isotropic, homogeneous, linear- 
elastic, continuous bulk medium with a smooth crack, at first glance it appears 
that there are still many variables which could infiuence the propagation of the 
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crack; for example, the load on the surface of the body, type of material, tem- 
perature, ambient atmosphere, and the stress and fracture history. However, 
we are concerned here with the propagation of a crack given the elastic fields 
in the body, not the precise conditions that produced those fields. Also, many 
variables, such as stress history and temperature, can be included implicitly in 
the materials constants, which we also assume are known. This means that for 
this work, the variables of relevance are the elastic fields near the crack tip, 
materials constants, and the current configuration of the crack. 

Now, it is well known that for cracks in linear elastic media, the stress field 
near the tip of the crack — which is the only area we expect to influence crack 
growth — obeys a ~ power law in the distance r from the crack tip in the 

plane. There are three modes of cracking (see figure la, b, and c); each has a 
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characteristic known angular dependence, which can be written 

^i.(r,^) = ^/5(^), (1) 

where is the angle made with n, the direction of cracking (see figure Id), 
and the three are the mode I, II, and III Stress Intensity Factors (SIFs), 
numbers characterizing the strength of the stress singularity for each cracking 
mode. There are also similar expressions for the displacement field near the crack 
(see reference 17). From equation 1, we can see that the relevant information 
from the elastic field for a two dimensional crack propagation problem, instead 
of being given by a displacement vector at every point in the plane, is reduced 
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to just three numbers, the SIFs. In weakly three dimensional problems, each 
non-equivalent plane has three SIFs which characterize the stress, so that the 
relevant information from the elastic field is given by three numbers at each 
point on the crack front, instead of a displacement vector at every point in the 
body. 

Now, the configuration of the crack could be specified by giving equations 
for the two surfaces of the crack, but it can also be reduced to a smaller amount 
of information. This comes about because part of the information on the crack 
geometry is contained in the elastic fields, and the part of the crack surface far 
from the crack tip cannot influence the growth of the crack beyond affecting the 
elastic fields. So, the most relevant information about the crack geometry, aside 
from the elastic field, is given by the crack front curve, x{X), and the vector 
n(A) giving the current direction of crack growth, where A parameterizes the 
crack front curve. In figure Id we show the three unit vectors associated with 
this description: t{X) = the tangent to the curve, with s the arc length; 
n(A), the direction of crack growth, perpendicular to t{X) (so that t{X) and 
n(A) define the crack plane at point A), and b{X) = i{X) x n(A), the normal 
to the local crack plane. Note that since the material containing the crack is 
isotropic, the coordinate system defined by these unit vectors is the only one 
physically relevant to crack growth, and all other quantities (such as the SIFs) 
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(primary stress is b-n). (c) A Mode III crack (primary stress is 
b-t). The local stress field around a crack can always be decom- 
posed into a linear sum of the three modes; the decomposition 
can be found from the modes' different symmetry properties (see 
appendix B). (d) The vectors associated with a point on the crack 
front: t is the tangent to the crack front curve; n, perpendicular to 
t and in the crack plane, is the direction of crack growth; b = txn 
is the normal to the crack plane. 



9 

are understood to be defined in this coordinate system (see figure 1). Also, in 
two dimensional cracks, t, h, and b are all constant along the crack front. 



IV. The crack growth law 

Now we are ready to derive a crack growth law in our relevant variables: 
materials constants, the SIFs Ka{X), and the unit vectors n(A), i{X), and b{X). 
That is, we are ready to derive an expression for the time evolution of the crack 
front curve x(A) as an expansion in ^ = (||) ^ the gradient along the 
crack curve, of the relevant variables. (In two dimensions, ^ = 0.) Noting that 
the crack surface is smooth for all time, we know that the time derivative of x 
must lie in the crack plane, so we can write 

= ^(A, t)n{X, t) + w{X, t)t{X, t), (2) 

where v is the crack velocity, and ta is a non-physical function which can be 
chosen freely to determine how A parameterizes the physical crack surface. That 
is, a particular choice of w determines a gauge for A (see appendix A); in two 
dimensions, we take w = 0. 

To determine the growth of the physical crack, then, we need to find ^(A, t) 
and n(A, t) in equation 2, in terms of the relevant variables. Noting that ?7.(A, 0) 
is given, this means that to find n(A, t), we need to find the time derivative 
Now, since fi is a unit vector, its time derivative must be perpendicular to itself; 
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also, by definition n is perpendicular to £ = This gives us a constraint on 
the equation of motion for n, obtained by setting ^ (n • f ) = : 

On ^ dv di ^ , , 



Note that in two dimensions the right hand side of this equation vanishes, since 

Another constraint on and one on the crack velocity v are obtained 
from symmetry. We consider symmetry operations, centered at some point on 
the crack front, which leave the unit vectors at that point fixed and refiect or 
rotate the material, preserving the physical properties that t is the tangent to 
the crack front curve and n is the direction of cracking at that point. (This 
is equivalent to leaving the material fixed and transforming the coordinates.) 
There are two such independent operations: (a) 180° rotation about n, and 
(b) reflection in the n-t plane (the crack plane) (see figure 2). Under both of 
these operations, the physical law for v must remain unchanged; that for^-b 
must change sign under both operations (see figure 2). 

Let us now examine the case of a crack in two dimensions, where the only 
non-vector quantities we can form are combinations of the SIFs and material 
constants. Under symmetry operation (a), the material constants, Kj, and Km 
remain the same, while Ku changes sign; under operation (b), both Ku and 
Km change sign, and everything else remains the same (see appendix B). Since 
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Fig. 2 Transformations (a) rotation about n and (b) reflection in the n-t 
plane on (c) the untransformed crack. To make a physical crack 
growth law, the crack velocity must be unchanged in both cases 
(if it changed sign, the crack would heal); ^"'gl'^^ ■ b must change 
sign in both cases, since the direction of b is unaffected by the 
transformation, while the physical normal to the crack plane shifts 
to the opposite direction. 
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the three SIFs have different transformation properties, and only Ku transforms 



where both v and / are functions of materials constants, Ki, K^^, and i^nj. (The 
dependence of / and v on Ku and Kui must be quadratic to insure invariance 
under both symmetry operations.) The minus sign makes / > correspond to 
the observed direction of crack growth under mode II loading; this is discussed 
in the next section. 

In three dimensions, the gradient ^ is not strictly zero, and there are 
therefore non-vector quantities besides the SIFs and materials constants which 
can be formed from the relevant variables. Up to first order in these are 
listed, along with their transformation properties, in table 1 of appendix B. 
From the transformation properties of these quantities, and from the discussion 
in appendix A on gauge invariance, we can see that the physical crack growth 
law, to first order in has the general three-dimensional form: 



hkef-S, 



we see that the two-dimensional crack growth law must have the form: 




dx 



(4) 



dx 
'di 
dfi 
'di 



vn + wi 




ds 



(5) 
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where the /, g^, and hij are functions of materials constants, Kj, Kij, and 
K^ii, and w is a function of materials constants, Ki, K^^, -f^ni? (-f^ii-^iii^^)) 
(Km^), (Ku^), {Ku§ ■ b), (|f • £), and {Kui^ ■ b), since these are all 
the quantities up to first order in ^ which are invariant under both symmetry 
operations. 

V. The undetermined functions in the crack growth laws 

Now, in equations 4 and 5 we have general forms for the crack growth law 
for two and three dimensional cracks. These equations contain many unspec- 
ified functions, which must be determined from considerations other than the 
symmetry considerations we used to find the general forms. Perhaps the most 
important of these functions is the crack growth velocity, v, a function of ma- 
terials constants, iCi, Kj^j, and K^^^ in two dimensions. This function has been 

1 3 

measured for mode I cracks and usually has the form shown in figure 3a, with 
zero velocity for Ki below some value Kj^, which depends on the material, and a 
monotonically increasing velocity above Ki^. This schematic form for the veloc- 
ity has also been found in a theoretical calculation for a viscoelastic system. ■'^^ 

1 S— 9 1 1 Q— 9 1 

However, both in theory and in very clean experimental systems, the 

crack velocity has the form in figure 3b, with a negative velocity (crack healing) 

when Ki < Ki^. This means that the crack velocity is a continuous function 

which passes through zero at Kj^, so that for SIFs near Kj^, (i. e. small crack 
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velocities), we can expand v as v{Ki) pa vq ^'^^^'^ , with vq a material dependent 
constant. For modes II and III cracks, as well as mixed mode cracks, since the 
elastic energy released per unit area of crack surface (the "energy release rate" ) 
is proportional to {Kf + Kf^ + ^j-z^), where p is Poisson's ratio, we expect 
that a crack velocity function valid for all modes of cracking can be expanded 



as: 



v{Ki, Kii, Km) pa vq — , (6) 

Kc 

where K = {K^ ^K^i + ^-l) ) ^ ■ Also, note that for fatigue cracking, where our 
growth laws must still hold (on time scales long compared to the load cycle), 
the crack velocity generally does not go to zero sharply at Kc-, but has a more 
gradual turn-on behavior (see figure 3a). 

Now, we saw in the previous section that the crack velocity in three dimen- 
sions can also depend on gradient quantities, besides the SIFs; the dependence 
of the crack velocity on these quantities has not been measured, to our knowl- 
edge. However, there is no reason to suppose that the dependence on these 
quantities has special behavior (e. g. zero-crossing or very strong dependence on 
SIFs) near K^i the value of the SIF where the crack velocity becomes positive. 
So, for small velocities, where K Kc-, we can approximate the dependence of 
the crack velocity on these quantities by a constant function (i. e. its value at 
Kc). 



V 
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Fig. 3 (a) Form of the c:-ack velocity v as a function of stress intensity- 
factor K in ordinary crack experiments, where the crack velocity 
is zero below a cri leal value Kc of the stress intensity factor, and 
then has a sharp turn-on. The dotted line shows the behavior 
under fatigue, where the crack velocity increases more gradually, 
(b) Form of the crack velocity in very clean experiments, where 
crack healing can take place. 
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Similarly, we expect that the seven functions f,ga, and hij in equations 4 
and 5, which are allowed by symmetry to be functions of materials constants, 
Ki, Kj^j, and i^ni, can be approximated as constants when the velocity is small. 
This means that to find the material-specific form of the crack growth law, for 
small velocities, it is a reasonable approximation to measure only the linearized 
dependence of the crack velocity on K, and the constant parts of /, Qa, and hij. 



VI. Predictions of the crack growth laws 

Let us now examine the two dimensional crack growth law, equation 4: 

dx 

(4) 

When Kii = 0, this equation says that the crack grows in a straight line (since 
^ = 0), in agreement with the "principle of local symmetry" generally used 
to predict crack growth in two dimensions. However, the principle of local 
symmetry also says that Ku = is maintained at all times by the propagating 
crack — in effect, that the crack curves in such a way as to keep Ku = 0. Our 
law, in contrast, says that it is only a non-zero Ku which can make the crack 
curve, but that (with / > 0) the crack curves in such a way as to make Ku 
smaller (see figure 4). 

Now, we can resolve the differences between the principle of local symme- 
try and our crack propagation law by writing the crack velocity v as the time 
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(a) 



(b) 



Fig. 4 Schematic of observed crack j growth in mode II, where the crack 
curves so as to reduce the niode II stress, leavmg only mode I 
stress, (a) Our picture, where the crack curves gradually to the 
direction where Ku = 0, on a length scale of j^. (b) The tradi- 
tional picture, where there is a sharp kink to the direction where 
Kii = 0. Note that in the / — > oo limit, the two pictures agree. 
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derivative of the crack length ^ , writing ^ = §f ^ : and by writing n and b in 
terms of the angle 6 that n makes with the x-axis. With these changes, equation 
4 becomes: 

dx 

— = (cos6,sm9) 
da 

I 

In principle, /, v, and Ku are functions of x and 9. However, in the case of 
a small amount of growth at the end of a long crack, we expect / and v to 
be nearly constant as the crack grows, since the SIFs only change by a small 
amount during the growth (see appendix C.) Also, when 6 differs from the angle 
that makes Ku = by only a small amount A9{x), we can approximate Ku as 

A9{x) 



Ku{x) = K,{0)^[l + o{l)], (8) 



where a is the length of the original long crack (see appendix C). Using equation 
8 in equation 7, we see that 



da \ 2v 

which we can immediately solve, taking / and v constant, to find that 



/KT(0)a 

A9{a) = e —. (10) 



That is, if we start a crack with a small deviation from the direction predicted 
by the principle of local symmetry, then our crack propagation law says that the 
deviation decays with a characteristic distance of j^. This length scale must 
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be very small, about the size of the non- linear process zone or the smoothing 

length, because these microscopic lengths are the only length scales that appear 

in two-dimensional crack problems (see section II). In the limit that the length 

scale is zero, or / — > cxo, we can now see that our crack propagation law for two 

dimensions agrees with the principle of local symmetry. 

Now, let us move to consideration of equation 5, 
dx 



„ -vn + wt 
at 



dn 
'dt 



dv dt ^ 
— +w— ■ n 
ds ds 



t+ 



(5) 



-/All + giKiii— h guKiiKiu— h gm^ h 

ds ds ds 

htb-^ • b + hntKii— ■ t + {hnbKiiKiii + w) — -b 
ds ds ds 

our crack propagation law for three dimensions. First, note that the principle of 
local symmetry term, —fKu, appears in this law, just as in two dimensions, and 
we do not expect / to be different here from its value in two dimensions (that is, 

is the inverse of a microscopic length). In contrast to /, the other functions 
appearing in equation 5 do not contain any length scales — both and 
are dimensionless. The length scales in these terms come from the gradient 

if the dimensionless forms of g and h are of order 1, these terms act over 
the length scales of the gradients. As noted in the previous section, we do not 
know ga or hij from symmetry principles; simple experiments to measure the 
sign and magnitude of these material-dependent functions, and their physical 
interpretations, will be discussed in a separate paper. 
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Our crack growth law, since it contains terms besides the principle of local 
symmetry term, predicts three-dimensional behavior beyond the scope of the 
principle of local symmetry; one example of this is the so-called "factory roof 
structure seen in mode III cracks. We will explore this and other predictions 
of our crack propagation law more fully in a separate paper, and also discuss 
the stability of straight cracks to wavy perturbations. 

VII. Conclusions 

We have seen that from symmetry principles, we can derive a crack growth 
law for both two and three dimensional geometries which agrees with the prin- 
ciple of local symmetry^^ in the limit that the microscopic length scales in the 
crack problem are truly zero. Our law also predicts behavior seen in three di- 
mensions that is not predicted by the principle of local symmetry; more detailed 
predictions of our crack growth law in three dimensions will be explored in a sep- 
arate paper. The laws we derived, equations 4 and 5, contain several functions, 
such as the crack velocity, which are not determined by symmetry alone and 
must be measured in controlled cracking experiments or atomic simulations; we 
will discuss experiments designed to measure the sign and order of magnitude 
of these functions, and their physical interpretations, in a third paper. 

By using symmetry principles, separation of length scales, gauge invariance, 
and gradient expansions, then, we have derived effective, macroscopic equations 
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governing the growth of cracks in three dimensions. We have course-grained 
the problem so that microscopic details — such as atomic bond breaking, crys- 
talline grain morphology, deformation near the crack tip in response to strain, 
and surface effects — are on such small length scales that they cannot affect the 
macroscopic crack growth. Understanding the microscopic origins of our effec- 
tive growth equations, and describing crack growth on very small length scales, 
where our crack growth law is not valid, will demand calculations that include 
these microscopic details. 

Acknowledgments 

We acknowledge the support of (US) DOE Grant #DE-FG02-88-ER45364; 
JH was funded by an NSF Graduate Fellowship and a DOEd National Needs 
in Materials Physics fellowship. We would like to thank P. Wash Wawrzynek, 
Dave Potyondy, and Tony Ingraffea for introducing us to the problem of crack 
propagation, for many discussions about the subject, and for allowing us to use 
their two-dimensional finite element analysis code (FRANC); Jim Krumhansl 
for his continuing support and critical evaluation of this research; and Andy 
Ruina, Robb Thomson, and Jim Rice for their comments on earlier versions of 
this work. JPS would like to thank Stephen Langer for introducing him to the 
gauge invariance of curves. We would also like to thank the Technical University 
of Denmark and NORDITA for support and hospitality. 



22 

Appendix A: Gauge symmetry and cracks 

There are many cases where the natural mathematical description of a prob- 
lem introduces fictitious degrees of freedom with no physical relevance. The 
most well-known example is in electromagnetism, where all physical quantities 
are unchanged when the gradient of an arbitrary function x(r) is added to the 
vector potential A. The transformation A + Vx is an example of a "gauge 
transformation;" the invariance of physical quantities under such a transforma- 
tion is called "gauge invariance." The strong and weak forces of particle physics 
also have gauge transformations associated with them.^^ In general relativity, 
the choice of coordinate system for space-time is arbitrary; this gauge invariance 

or 

can be used to derive momentum and energy conservation.^ Another use of 
gauge symmetry is in site-disorder spin glasses, where a gauge transforma- 
tion is used to show that the certain forms of disorder do not result in spin glass 
properties, but in ferromagnetic behavior. 

The term "gauge" is particularly appropriate for the gauge symmetry of 
our problem, ' where the parameterization A of the crack front curve x{X,t) 
is arbitrary: how one "gauges" (measures) the points along the curve cannot 
affect the growth of the crack. There are two different types of gauge symmetry 
for cracks. The first type is the freedom to change the parameterization at any 
one time, which we call the "one-time gauge symmetry." The second is the 
freedom to choose how the parameterization at some time is related, through 
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the growth equation, to the parameterization at a later time; we call this the 
"time-dependent gauge symmetry." Crack growth laws must satisfy both gauge 
symmetries; that is, neither the one-time nor the time-dependent gauge trans- 
formation can change the physical crack growth equation. 

A crack growth equation that satisfies the one-time gauge symmetry must 
not have any direct dependence on the value of A at a point on the crack front 
curve, but only depend on physical quantities evaluated at that value of A. 
Also, derivatives along the crack front cannot enter the growth equation as 
but must instead be in terms of the arc length s, because ^ is gauge invariant. 
Since we have written our crack growth law in terms of ^ , and have not included 
explicit A-dependence, it does satisfy the one-time gauge symmetry. 

The time-dependent gauge symmetry is slightly more complicated. An 
equation for ^^^'^^ , which is how we have chosen to write the crack growth 
equation, tells us how a point with parameter value A evolves in time. This 
means that the time evolution of the parameterization is implicit in our for- 
mulation, and the time-dependent gauge transformation changes the growth 
equation (unlike the one-time gauge transformation). This change happens in a 
well-defined way: if we have some crack growth law in terms of a parameter A, 



dx{\, t) 

dt 
dn(X,t) 

dt 



= Ci+Db, 



= Ah + Bi 
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where the right hand sides are imphcit functions of A, when we introduce a 
time-dependent gauge transformation to a new parameter fj,{X, t), then the crack 
growth law becomes 

at 9/1 at 

Ot OfJ, ot 

with the right hand sides now imphcit functions of /j,. Writing ^ = ^ (with 
s the arc length), defining a new function w = and using the definition 

of £ = 4f , we can write this as 



ds ■ 

dx{fj,, t) 
dt 



An+{B + w)i 



?i^=Ci+Dh + J^ (13) 
ot OS 

= {C + w^-t)t+{D + w^-b)b. 

OS OS 

There are three particular time-dependent gauges that we have found to be 

of special use in our study of crack growth. First, there is the "reference gauge," 

where B + w in equation 13 is zero. In this case, curves of constant parameter 

value A„ are the integral curves of n, and the growth equation can be written: 

dx{Xn,t) 

(14) 

dn{Xn,t) _ _dVf ^ 

dt ~ as ^ 

Here the term comes from the requirement that n-i = Ohe preserved at all 
times; the function D is free, as far as gauge symmetry is concerned. This is the 
only possible form of the growth law when ^ is along n, and we use this as our 
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reference gauge for discussing other time-dependent gauges. In fact, comparing 
equation 14 to equations 11 and 13, we see that the growth law in a general 
time-dependent gauge can be written as: 

(15) 

with w = as above. Note that the function w characterizes the time- 

dependent gauge used in the crack growth law, while D and v are the physical 
functions describing the crack growth. 

Another particular time-dependent gauge of interest is the "arc length 
gauge," where the parameter A is always equal to the arc length s along the 
crack front curve, measured from some starting point (such as the edge of the 
body). Since the arc length along the crack front is s{X) — J^^ \ \^\\dX, if we 
begin in the arc length gauge (by making a one-time gauge transformation), 
then we can remain in that gauge by choosing 

^A,) = j\{X')^.i{X')dX' (16) 

in equation 15. Although arc length is physically the most natural parameteri- 
zation for a curve, the arc length gauge is not usually very convenient, as w{Xs) 
is a non-local function — as the crack grows, if the arc length of a section near 
some point Ai stretches (or shrinks), A^ must shift upwards (downwards) for all 
points with A^ > Ai. 
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A third time-dependent gauge, the "z gauge," is useful for cracks which 
have crack fronts which always point nearly along some axis, which we take to 
be the 2;-axis. In this case, it is natural to use a gauge where the parameter 
A of the crack front is the ^-coordinate (as long as the crack front ^(A^) is a 
single- valued function). To achieve this, the crack growth equation for x must 
have zero ^-component; that is, we must choose 

w{X,) = -v{X,)'^. (17) 

Similar ideas can be used to make special gauges for cracks which are nearly 
circular, parameterizing with the angle 9 from the a;-axis, for instance, and for 
other common crack geometries. 

Appendix B: Symmetry operations on the relevant variables 

In this appendix, we examine the two symmetry operations (a) 180° rota- 
tion of the material about n at some point Aq, while keeping the coordinates 
fixed, and (b) reflection of the material in the n-t plane (the crack plane) at 
Aq, while keeping the coordinates fixed (see figure 2), as applied to the relevant 
variables for crack growth and their derivatives. Note that only the signs, and 
not the magnitudes, of the variables and derivatives can change under these two 
operations. Also note that the signs of the SIFs are defined in terms of the three 
unit vectors, and that the gradient operator ^ along the crack front curve goes 
in the direction of t. 
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Now, Ki transforms like the Ufeb stress component, Kn like anh-, and Km 
like (Tfot (see figure 1). Symmetry operation (a) takes material at 3:5, a^^), in 
terms of the coordinates with axes n, 6, and t and origin at Aq, to (a;„, — X5, — a;^); 
symmetry operation (b) takes material at x^, Xf) to (x^, — X5, Xt). Therefore, 
— > under both (a) and (b), Kn — > — -fCn under both, and Km — > Km 
under (a) and —Km under (b). 

Now let us consider a case where one of the SIFs, before transformation, is 
greater, in absolute value, for > than for xt < 0, so that > 0. Then, 
under transformation (a), the material at > with the greater \K\ moves 
to < 0, so that < 0; transformation (b) leaves xt unchanged, so that 

> 0. This can be combined with the transformation properties of the SIFs 
themselves to give the transformation properties of the gradients of the SIFs; 
the results are in table 1. 

We also need to consider the transformation properties of gradients of the 
unit vectors — quantities of the form ^ • b, where a and b belong to {n,b,t}. 
Noting that 

d(a-b) ^ da i db ^ 

for all a and &, since n, 6, and t are mutually orthogonal unit vectors, we can 
see that there are only three independent quantities to consider, which we take 
to be ^ • £, • ^, and || • b. From figure 5, we can see that the transformation 
properties of these three quantities are as shown in table 1. 
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Quantity 


loU Tt rotation 


(b) 71 — t reflection 


Ki 


1 

+ 


1 

+ 


An 








1 




dKi 
ds 




1 

+ 


dKu 

ds 


1 

+ 




dKiu 
ds 






dh f 
ds ' 


+ 


+ 


ds 


+ 




ds " 







Table 1 Transformation properties of relevant variables and their deriva- 
tives under two symmetry operations; see figure 5 and text. Note 
that products of these quantities transform as the product of the 
transformation properties (e. g. {KuKm) is — under (a) and + 
under (b)). 

Appendix C: Approximation of Ku in two dimensions 

In this appendix, we derive equation 8, 

M{x) 



Ku{x) = Ki{0)- 



l + O 



(8) 



which gives Ku after the crack has grown a distance x from the end of a long 
crack of length a, in terms of the deviation A9{x) of 9 (the angle that n makes 
with the X-axis) from the angle that makes Ku = 0. First, we can use the 
results of Cotterell and Rice^^ to find that as a function of the x-coordinate of 
X, measured from the end of the original long crack. 



Ku{x) = Ku{0) + le{x)K,{0) -(IYt r . 



-dx', (19) 
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Fig. 5 (a) Top view of tlie surface of a planar crack witli a curved crack 
front, which gives a non-zero value for ^ • t, with both ^ • b 

and ^ ■ b zero (since both h and i are always in the same plane, 
perpendicular to b). Under rotation of the crack about n at a point, 
keeping the vectors fixed (symmetry operation a), the value of ^-i 
stays the same; under reflection of the crack in the n-t plane at a 
point, again keeping the vectors fixed (symmetry operation b), it is 
also invariant, (b) A non-planar crack with non-zero ^ • b, where 
both ^ ■ i and ^ ■ b are zero (since b and h are always in the same 
plane, perpendicular to t). Under symmetry operation a, ^ ■ b 
is invariant; it changes sign under symmetry operation b. (c) A 
non-planar crack with non-zero ^ • b, where both ^ • i and ^ • b 
are zero (since b and t are always in the same plane, perpendicular 
to 77.). If • & changes sign under symmetry operations a and b. 
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where T is the non-singular tensile stress at the end of the crack. So, when 
the principle of local symmetry is satisfied, 6{x) has the value which makes 
Kii{x) = 0; if ^ differs from this value by a small amount A^(a;), and if we take 
T — ; with the appropriate geometrical factor 6, then we find that 



Kii{x) = Ki{0) 



Ae{x) 



I -Ah 



2xY r A9{x'] 



Tra 



AO{x) \x — x' 



X 



2 dx' 
~2x 



(20) 



Now, if ^^Q^J^ were constant, then the integral on the right hand side would be 
equal to 1; small variations of AQ{x) from a constant function leave the integral 
approximately 1. Noting that the integral is multiplied by (f ) ^ , which is small 
by assumption, we can therefore approximate K\\ as 



K^iix) = Ki{0) 



A9{x) 



(8) 



which is equation 8. We also note here that under the approximations of this 
appendix, and with the results of Cotterell and Rice,''^''^ Ki is approximately 



constant. 
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